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Functions of Matrices. 

By H. B. Phillips. 



1. It is the purpose of the present paper to study the functions repre- 
sented by polynomials or convergent series in a matrix or a finite number of 
matrices. As the work is concerned mainly with the roots of the matrices, the 
fundamental facts about the roots are first briefly developed.* 

By a matrix of the n-th. order is meant a square array of n 2 elements a ik , 
i, k = l, 1, . . . ., n, 

a n a 12 . . . . a ln 
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a nl a n2- • • - a n 

The matrices considered in this paper will all be of the same order. 
The determinant 

a n a 12 . . . . a ln 

$21 $22 • • • • $2» 
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= a, 
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is called the determinant of A. When this determinant is zero, A is called 
singular. 

The sum of two matrices J. = ||« d || and i?=||&<,,|| is the matrix 

A+B=\\a ik +b ik \\. 

More generally, if a, and (i are numbers, hA+ [iB =\\ha ik -{- [ib ik \\. A matrix is 
zero when and only when all its elements are zero. 

The product AB of -4=||o rt || and B=\\b ik \\ is the matrix 

n 

AB=\\a ile \\, where a ik =Ha ii b jk . 

* The general theory of matrices is given in Bocher's " Introduction to Higher Algebra." A very 
complete bibliography of literature before 1907 is given in James Byrnie Shaw's " Synopsis of Linear 
Associative Algebra," published by the Carnegie Institution of Washington. 
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The products AB and BA are not in general equal. If these products are 
equal, A and B are called commutative. 

The products of three or more matrices are associative, that is, 

(AB)G=A(BC)=ABC. 

The determinant of a product of matrices is equal to the product of their 

determinants, for example, | ABC | = | A | • \B\ • \C\. 

The matrix 

1 0....0 

1 0....0 

0....1 

in which the elements a u are all unity and the others zero, is called the unit 
matrix. It is easy to see that AI—IA = A. 

If the determinant of A is not zero, there is a matrix A" 1 called the 
reciprocal of A, such that AA~ 1 —A~ 1 A=I. In this case, if AB—CD, 
B^A^CD, i. e., we can divide by a non-singular matrix. 

It is to be noticed that in both multiplication and division the operation 
can be performed on the right or on the left. B is multiplied on the left by 
A if the result is AB, and on the right by A if the result is BA. Similarly, B 
is divided on the left by A if the result is A~*B and on the right if the result 
is BA- 1 . 

2. Identical Equations. — The matrices A, B, . . . ., P satisfy the equation 

A + B+ +P = 

if their sum is a matrix \\a ik \\ with elements u ik all zero. Cayley first observed 
that a matrix of the n-th order satisfies an algebraic equation of the n-th 
degree. This may be considered a consequence of the following theorem : 

Theobem I. Let A = \\a ik \\, B — \\b ik \\, . . . ., P=\\p ik \\ be a finite number of 
matrices of the n-th' order and let u ik =^a ik +(ib ik + . . . . -{-opa, X, (i, ....,p 
being numerical parameters. If A', B', . . . . , P' are commutative with each 
other and satisfy the equation 

AA' + BB' +.... +PP'=0, (1) 

they will also satisfy the n-th degree equation 

\a {k A' + b ik B'+....+p ik P'\=0, (2) 

obtained by replacing %,, (i, . . . ., p in |a iA .| =0 by the matrices A', B', ... ., P', 
respectively. 
34 
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This theorem enables us to replace a homogeneous linear equation with 
matrix coefficients by one of the n-th degree with numerical coefficients. The 
latter equation will usually have solutions, however, that do not satisfy the 
former. 

To prove the theorem, let E ik be the matrix of the n-ih order with all its 

elements zero except that in the i-th row and k-th column, which is unity. It 

is readily seen that 

E {j E jk =E ik , E {j E lk =0, Ixpj. (3) 

The matrix A can be written „ 

A = S a ik E ik , 

i,h=l 

and so (1) is equivalent to 

2 E ik (a ik A' + b ik B'+ .... +p ik P') =0. 

Multiplying this equation on the left by E n , E n , . . . ., E ln , respectively, and 
using (3), we get 

E n (a 11 A'+ .... + Pll P') +E n (a n A'+ .... +p 12 P') + .... +E ln (a ln A'+. . . . +p ln P')=0, 

E u (a n A'+ .... +PnP') +E 12 (a 2i A'+ ••■■ +fe-P') + • • • • +E ln (a in A'+ .... +p 2n P')=0, 

• • > 

E 11 (a nl A'+ ■■•• +P n iP') +E n (a n2 A'+. . . . +p n JP') + .... +E ln (a nn A'+ .... +p nn P')=0. 

Since A', B', . . . . , P' are commutative, we can eliminate E 12 , E w , . . . . , E ln 
by multiplying these equations (on the right) by the cof actors of the corre- 
sponding elements in the first column of the determinant 

\a {k A'+b ik B'+....+p ik P'\ (4) 

and adding. The result is 

E u \a ik A' + b ik B'+ .... +p ik P'\ =0. (5) 

If we expand (4) and combine terms, the result is a matrix. Equation 
(5) shows that all the elements in the first row of that matrix are zero. Tor, 
the product of E u and any matrix ||a £ J| has as first row a n , a 12 , . . . ., a lB . 
If then ^nlla^H is zero, the elements a u ,a 12 , ....,a la must all be zero. 
Similarly, we could eliminate E n , E 18 , . . . . , E ln and so get 

E 12 \a ik A'+b ih B'+ .... +p ik P'\ =0, 

showing that all the elements in the second row of the matrix represented by 
(4) are zero. By continuation of this argument we conclude finally that 

\a ik A'+b ik B'+....+p ik P'\=0. 
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3. Characteristic Equation.^Let I be the unit matrix. Applied to the 
identity AI — 7*4=0, Theorem I gives 

ai\I — A, a 12 x, .... a ln I 
02i~i , a^l — A. , 



. . . a 2 J 



a n J- 



=0. 



This is an equation of the w-th degree in A called the characteristic equation 
of A. Arranged in descending powers of A (with a change of sign if n is 
odd), it takes the form 

<t>(A)=A»+a 1 A— *+.... +aj=0. (6) 

Let % , a 2 , . . . . , a n be the roots of the equation 

$(r) =r n +aif ffl-1 + + a„=0. 

Then $(r) can be factored in the form 

4> ( r ) = (r— o x ) (r— a 2 ) .... (r— a n ) . 

Since this is an identity in r, the coefficients of each power of r on the two 
sides of the equation are equal. It will then still hold when r is replaced by 
A, and a,, . . . ., a„ by aj, . . . ., aj. Hence 

${A) = (A—aJ) (A—aJ) (A—aJ). 

The numbers a 1} a 2 , . . . ., a n are called roots of the matrix A. These 
roots satisfy the equation 



a n - 

«21 



r a 12 



a. 



,—r. 



"In 



*nl 



This expresses that the determinant \A— rl\ is zero. The roots of A are 
therefore the values of r for which the determinant \A — rl\ is zero. 

A matrix of the n-th. order may satisfy an equation of lower than the n-th 
degree. The equation of lowest degree satisfied by a given matrix will, how- 
ever, be unique. For, if A satisfies two equations of the m-th degree, we can 
eliminate A m from them and so find an equation of lower degree satisfied 

by A. Let 

+ (A) =A m +c 1 A m ~ 1 + .... +cj=0 

be the equation of lowest degree satisfied by A. This is sometimes called the 
reduced equation for A. It is clear that 

^(r) =r m +e 1 r m ~ 1 + . . . . + c m 
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is a factor of <p(r). For, by division, we get 

cp(r)=Q(r)4,(r)+R(r), 

where Q(r) is the quotient, and R(r), of degree less than m, is the remainder. 
Since this is an identity in r, we can replace r by A and so get 

^(A)=Q(A)^(A)+R(A). 

But ty{A) and ^{A) are both zero. If then R(r) were not identically zero, 
R(A) =0 would be an equation of lower than the m-th degree satisfied by A. 
Hence R(r)=0 and $(r) = Q(r)^{r). 

This shows that all the roots of the reduced equation are roots of the 
characteristic equation. Conversely, all the roots of the characteristic equa- 
tion satisfy the reduced equation. For, since ^(r) is a polynomial, we can 
factor ^(r) — ^(s) m the form 

4>(r)—4'(s) = (r—s)P {r,s), 

where P(r,s) is a polynomial in r and s. Since this is an identity we can 
replace s by A and r by rl. Then, since 4 (A) = 0, 

4>(r)I=(rI— A)P(rI,A). 

Equating determinants of the two sides, we get 

[^(r)]"=|rZ— A\ • \P(rI,A)\. 

If now r is a root of A, \rl — A\ —0, and so 4 , (») =0 which was'to be proved. 
4. Associated Roots. — Theoeem II. If A, B, . . . ., P are commutative 
matrices with roots %, a 2 , ....,«„; b 1 ,b 2 , . . . .,b n , etc., those roots can be 
arranged in sets a if b i} . . . ., p i} i — 1, 2, . . . ., n such that, if f(a, b, . . . ., p) 
is any polynomial in a, b, . . . ., p, the roots of the matrix f(A, B, . . . ., P) are 

/(a,-, b if , Pi), i = l, 2, ,n. 

This theorem is due to Frobenius.* The following proof is somewhat 
more direct than the one given by him. Let f x (A, B, . . . ., P), f 2 (A, B, . . . . , P), 
etc., be polynomials and 

f(A,B, ....,P) = £\f i (A,B, .. .,P), (7) 

Xj, a. 2 > • • • -, ^k being arbitrary numbers. Since A, B, . . . ., P have character- 
istic equations of the m-th degree, there are only a finite number of linearly 
independent polynomials in those matrices. We may consider f t (A, B, .... , P), 
t = l, 2, . . . . , h as those polynomials, and so by a proper choice of 1 1} \, . . . . ,\ 
make f(A, B, . . . ., P) equal to any given polynomial in A, B, . . . ., P. 

* Sitzungsberichte Berliner Akademie (1896), p. 602. 
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By Taylor's theorem, 

f(A,B, ....,P)-f(a i ,b j , ....,p k )I 

= (A- ai I)L+ (B-b j I)M+ . . . . + (P-p k I)R, (8) 

where L, M, . . . ., R are polynomials in A, B, . . . ., P, a j( bj, etc. Form the 

product 

U[f(A,B, ....,P)-f(a t ,b„ ....,p k )I] 

of all the differences (8) obtained by letting i, j, . . . .,k vary independently 
from 1 to n. This product is zero. For when the product is expanded in 
terms of the quantities (A — aj), (B — bjl), etc., each term of the result not 
containing (B — bjl), . . . ., (P — p k I), j, . . . ., k being definitely assigned, will 
contain all the factors (A — a^l), i=l, 2, . . . .,n. Hence each term contains 
one or more of the products 

(A— aj) (A—a 2 l) (A— aj), 

(B-b 1 I)(B-b 2 I)....(B-bJ), 



(P- Pl I)(P-pJ)....(P-pJ). 

These all vanish because the result of equating any one of them to zero is the 
characteristic equation of one of the matrices. Therefore 

U[f(A,B, ....,P)-f(a { ,bj, ....,p»)I]=0 (9) 

for all values of the parameters \ , \ , • • ■ • > \ • 

Equation (9) is a polynomial equation in f(A, B, , . . ,, P). If the equa- 
tion of least degree satisfied for all values of ^ , /l 2 > • • • • > \ hy that matrix is 

<M*)=o, (10) 

$(#) will be a divisor of 

II|>— f(a if bj, ,p k )]. 

The roots of (10), i.e., the roots of f(A,B, . . . .,P),are all therefore obtained 
by assigning proper values to i, j, . . . . , k in 

/(«<, bj, ,p k ). (11) 

We can therefore find n sets of constants a* , /^ , . . . . , p,- (namely, the sets of 
values a if bj, . . . ., p k in (11)) such that the roots of f(A, B, . . . ., P) are 

/( ai ,&, ...., 9i ), »=1, 2, ....,«. (12) 

In particular, if \, 2. 2 , . . . ., "k k are so chosen that f(A, B, . . . ., P) =A, we 
shall have /(a 4 , &, . . . ., p.) =a 4 . The roots of A are therefore a x , a 2 , . . . ., a n . 
Similarly, the roots of B are /? x , /? 2 , . . . . , /? n , etc. By a change of notation 
we can then make a ( =a 8 -, ^ i =b i , etc. The roots of f(A,B, . . . ., P) are then 
f( a i>bi, . . . ., p^, which was to be proved. 
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5. Canonical Form. — Two sets of associated roots o x , 6 1} . . . ., p x and 
a 2 ,b 2 , . . . ., p 2 will be considered different unless 

Let ai,bi, , p it i=l, 2, ,m 

be tbe different sets of roots of the commutative matrices A, B, . . . ., P. Let 
X, //, . . . . , p be numbers not satisfying any of the equations 

%a t +ltb t + +pjp<=to»+f*&»+ +pp k , *=£&. 

By Theorem II, the distinct roots of 7^A-\-^.B-\- . . . . -f-pP are 

Aa 1 +p& 4 + +pPi» i—1,2, , m. 

The factors of its characteristic equation are 

*A + tiB+ .... +oP- {Xa i + i ib i + • • • • +pPt)I 

^(A-aJ) +p(B->bJl+ +p(P- Pi I), 

and the equation of least degree satisfied for all values of %, p, . . . . , p is 

m 

n lUA-aJ) +(i(B-b { l) + .... +p(P-p,Z]'«=0, (13) 

the numbers r lf r 2 , . . . ., r m being the multiplicities of the roots. 

m 

Let * 1= n &(A-a t I) +^(£-6,1) + • • • • +p(P-p i I) ]'• 

«=2 

and, generally, let 4/ be the product of all the factors in (13) except 

[M^-«,D +!"(#-¥) + . • • • +p(P-p,I)]". 
It is clear that 

M»=0, »#:*, (14) 

for the product contains all the factors in the left member of (13). Let ^ be 
the function obtained by replacing X, fi t . . . . , p in ^ by W, /it', . . . . , p'. Con- 
sider the product 

[a'^-oj) + . . . . +p'(P_ Pl J) ]^. (15) 

When this is multiplied by ^[ the result is zero because 

WiA-aJL) + .... + p '(P-pJ) ]S4=0. 

Also the product of (15) and any one of the functions <4< 2 , ^S) . . . ., ^ is zero 
by (14). Hence 

[3/(4-0x1) + . . . . +p'(P-pJ)Y'Mtt J r4>2+ ■■■■ +4> m ) =0. (16) 
The roots of the matrix 

4-1+^2+ • . • . +<iv (17) 
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are obtained by replacing A, B, . . . . , P by a it b if . . . ., p t . When they are 
replaced by a lf 6 lf . . . ., p x the result is not zero because i^i is not zero, 
whereas 

4'2 = 4'3= ='4'm = 0- 

In a similar way it is seen that the other roots of (17) are not zero. We can 
therefore divide (16) by (17) and so get 

WU-<hI) + • • • . H-p'(P-PiI) ] r 4i=0. (18) 

Since X, fx.', . . . ., p' are arbitrary, if (18) is expanded in powers of those 
parameters, each coefficient will vanish. Hence 

(A— aJ)«(B— 6 X I)A. . . . (P-pJ)Hi=0 

for all positive integral values of a x , @ x , . . . . , p x such that 

ai+&+ +Pi=ri. 

In a similar way it is shown that 

(A-ajy^B-bjy*.. . . (P-pjyH><=o, (19) 

for all positive integral exponents such that 

a«+&+....+p,=r«. (20) 

Let *=*i+'k+....+<J'.. (21) 

As in case of (17) it is shown that none of the roots of (21) are zero. We 
can then define a set of functions $ { by the equations 

*,= ^, »=l,2,....,m. (22) 

It is clear that 

fc+fr+ .... +*.=/. (23) 

Also, by (14), 

M»=0, *=£&• (24) 

Multiplying (23) by ^ we therefore get 

*<=*<, *=1, 2, ...., w. (25) 

Let f(a,b, . . . ., p) be a polynomial in a, b, . . . ., p and 

2a+j8+....+p 

/../»,....,p(Oi 6i • ■■',?)= da a db" . . . .dp^ ^ a ' 6 ' • * • •' p) * 
By Taylor's theorem 

/(a, 6, ,p) =/(a„ 6 , p«) 

,yy y* /- i, , (a— a f ) a (g— ftV 

a p I a |p 
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Since this is an identity in the variables a, b, . . . . , p, we can replace them by 
A, B, , P. Then 



(26) 



f(A, . 




■,P)=f( 

+2.. 

a 


a iy 

p 


,Pi)l 

..,p\^i ) • • 


^ (A-ajy 
--,P<) ]JL " 


(P- 


-Piiy 


Let 




A, 


=^(A- 


-aji), ... 


^P^friP- 


-*>«*)• 






From (25) 


it 


follows 


that 
















Am-.- 


.Pf=fc(4-o,I)« 


l (P-W-. 


..(P- 


-*><*)'• 





Equations (19), (20), and (22) show that this is zero if a+/?+ .... +p>«V 

Finally, if we multiply (26) by &, sum for i = l, 2, , m, and use (23) we 

get 

/U, . . . ., P) = 2fc/(a„ . . . ., Pl ) + 2 2. . . .Ufa, M- ■ ■ -Pi) nr^Tn> ( 2 ?) 

the summation including powers A* . . . .P£ for which 

a + /3+....+p<r i . (28) 

Equation (27) gives a form in which any polynomial in the given com- 
mutative matrices A, B, . . . ., P can be expressed. We shall refer to it as the 
canonical form for a function of the matrices.* Its most important property 
is expressed in the following theorem : 

Theokem III. If A, B, ....,P are commutative matrices with corre- 
sponding roots a^bi, ....,p d and f(a,b, ....,p) is any polynomial, the 
matrix f(A,B, ....,P) can be expressed as a linear function of matrices 
depending only on A, B, . . . ., P, the coefficients in the linear functions being 
obtained by substituting each set of roots a^bi, . . . .,p i in f(a, b, . . . .,p) and 
in its partial derivatives of order lower than the multiplicity of the root 
a.a i +/K& i + .... -\-oPi in the equation of least degree satisfied for all values of 
\ ix, . . . ., p by KA+ixB+ .... -fpP. 

A case of particular interest is that of a polynomial such that all the 

coefficients f(a i} , p t ), /„,... .,/>(<*<> > Pi) in (27) vanish. Then, evidently, 

f(A, B, . . . ., P) =0. We can consider the values a it b { , . . . ., p t as coordinates 
of a point in hyperspace. In the case considered, the function f(a,b,....,p) 
has a zero of order r % at the point (a^bi, .:..,p t ). Therefore we have 
proved 

* The formula for a function of a single matrix with distinct roots was given by Sylvester, Comptes 
RendMS, Vol. XCIV (1882), p. 55. The case of a single matrix with repeated roots was given by 
A. Buchheim, Phil. M-ag., (5) 22 (1886), pp. 173-174. 
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Thboeem. IV. // the polynomial f(a,b, . . . ., p) has at each of the points 
(a it b iy . . . ., p^, i=l, 2, . . . ., m, a zero of order equal to or greater than the 
multiplicity of the corresponding root in the equation of least degree satisfied 
for all values of 1, [i, . . . ., p by %A+[iB-{- . . . . + pP then f{A,B, . . . .,P)=0. 

If the matrices fy, A^Bf. . . .P? in (27) are linearly independent, Theorem 
IV expresses the necessary and sufficient condition that f(A,B, ....,P) 
vanish. This is easily shown to be true for polynomials in a single matrix.* 
For polynomials in two or more matrices such may not be the case. For 
instance, A and B could be equal. Then A t and B { would be equal. 

6. Commutative Matrices not Expressible as Polynomials in the same 
Matrix.— The simplest illustration of commutative matrices is furnished by 
polynomials in a single matrix. If 

A = a 1 <pt>+0L2<P P ~ 1 + ■■■■ +«p Z > B = p i <p*+Pd>*- 1 + .... + &/, 
obviously A and B are commutative. If it were true that any two commuta- 
tive matrices could be so expressed,! by a repetition of the process, any finite 
number could be expressed as polynomials in the same matrix. The results of 
the preceding sections could then be more readily obtained by using these 
expressions. That such is not the case will now be shown by a simple 
example. Let 

A = 

By direct multiplication it is shown that 

A 2 =B 2 =AB=BA = 0. (29) 

Hence A and B are commutative. Suppose they are expressible as poly- 
nomials in q>. Since the characteristic equation for q> is of third degree, the 
expressions can be reduced to the form 

A = ai <p 2 +a 2 <p + aL 3 I, B^prf+prf+fal. (30) 

If a x and & are both zero, these polynomials are of first degree. If not, 
we form the expression 

(P^- ai P) 2 = [5 1 (a 2 4»+a 3 7) -a 1 (&$+ ft/)]*. (31) 

Equation (29) shows that this is zero. The expression 

[(3 1 (a 2 x+ a 8 ) — aiifax+ps) ] 2 

* See Metzler, Amehioax Joubnal, Vol. XIV (1892), p. 339. 

f Frobenius raised this question in the article to which we have previously referred, but stated that 
he had not decided whether it could be done or not. 

35 
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is not identically zero, however, because A and B do not satisfy an equation 

of the form ft L A — ^5 = 0. By the use of (31) we can then reduce (30) to the 

form 

A = a t $ + a 2 l, B = b 1 <p+b 2 I. (32) 

The coefficients a x and b x can not be zero because A and B are not multiples of 
the unit matrix. Hence, as before, we form the expression 

(b 1 A—a 1 B) 2 = (b 1 a 2 —a 1 b 2 ) 2 I = 0. 

This shows that b x a 2 — a t b 2 =0, and so, from (32), b 1 A—a 1 B=0. Since A is 
not a multiple of B this is impossible. Therefore A and B cannot be expressed 
as polynomials in $. 

7. Limits and Convergence. — A variable matrix ^> = ||r <s || is said to 
approach a matrix .4 = ||a 1Jfc || as limit if 

Lim r ik = a ik , i, k = l, 2, , n. 

In this case it is clear that each root of 4> will approach a root of A as 
limit. For the roots of q> and A satisfy the equations 

\rI-cp\=r»+p 1 r n - 1 + ....+p n =0 ) (33) 

\aI—A\=a n +a 1 a n ~ 1 + +a„=0. (34) 

The coefficients p x , p 2 > • • • • > P» are definite rational integral functions of the 
elements r ik of q>, and the coefficients a x , a 2 , . . . ., a n are the same functions of 
the elements of A. When <£> approaches A as limit 

Lim p<=a 4 , i — 1, 2, , n. 

Also, if the coefficient of the highest power (in this case unity) does not 
approach zero, the roots of a polynomial equation are continuous functions of 
its coefficients.* Hence as <p approaches A as limit the roots of (33) approach 
those of (34) as limits. 

Since different matrices can have the same roots, the roots of $ can 
approach those of A- when 4> does not approach A as limit. 

The sum q> 1 + $ 2 + ■ . . . +<p m of m matrices is equal to a single matrix 8 m 
If S m approaches a limit S when m increases indefinitely, the series 

fc+fc+.... +$.+ .... 
is said to converge arid have the sum S.\ In case of a multiple series, such as 

2fc.i,», h h k = \, 2, , co, 

* Weber's Algebra, Vol. I, § 44. 

fPeano treated the convergence of a matrix series by means of a modulus, Math. Annalen (1888), 
Vol. XXXII, pp. 450-456. See also G. A. Bliss, Annals of Mathematics (2) 6, pp. 49-68. 
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an order of summation is assumed to be included in the definition of the 
series. It is thus equivalent to a simple series. 

If A, B, . . . ., P are commutative matrices, by Theorem III, any poly- 
nomial f(A, B, ...., P) is expressible as a linear function of matrices 
depending only on A, B, . . . ., P, the coefficients being f(a { , b t , . . . ., p t ) and 
perhaps some of the derivatives fa, ft, . . . ., p(a if b if ....,Pi). If these 
coefficients approach definite limits when the number of terms in the poly- 
nomial increases indefinitely, the matrix f(A, B, . . . ., P) will then approach 
a definite limit. 

Theobem V. If A,B, . . . ., P are commutative matrices and f(a,b,....,p) 
an infinite series, the series f(A, B, . . . ., P) will converge and be represented 

by (27) if the series f(a { , b t , , p { ), f a ,p, „(«», b iy , p { ) in the right 

member of (27) converge. 

In this theorem it is understood that / aj3 „(a,b,....,p) is obtained 

by differentiating f(a,b, . . . ., p) term by term and arranging the results in 
the same order as the corresponding terms in f(a, b, . . . ., p). 

If f(a,b, . . . ., p) is an analytic function of the variables a, b, . . . ., p 
and A, B, . . . ., P an equal number of commutative matrices, we define 
f(A, B, . . . ., P) as the value (if it is definite) determined by (27). In order 
that the value be definite it is sufficient that f(a, b, . . . ., p) be analytic near 
each of the points (a t , b if . . . ., p t ) determined by the distinct sets of asso- 
ciated roots. For then the derivatives in (27) are also definite. If the 
function is many valued, a definite branch must be chosen at each of the points 
(aijbi, . . . ., p^. The same branch need not however be used for all. 

8. Taylor's Series. — Let f(z) be a function of the complex variable z 
analytic in the neighborhood of each root of a matrix A. If we choose a 
definite branch of /(e) at each root (but not necessarily the same branch for 
different roots) the functions f(A), f'(A), f"{A), etc. are determined by (27). 
Let Z be a matrix commutative with A and consider the series 

■HA,Z)=f(A)+f(A){Z-A) + ....+rU) K | m ^ + (35) 

Let s t and a t be corresponding roots of Z and A. The series (35) will con- 
verge if each root z t lies within the circle of convergence of f(z) with center 
at a t . For then the series 

/(«<) +/'(<*«) (e*-a<) + ■■■■ +f m (ai) ^T^" +••••, (36) 
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and the series obtained by substituting a { and s f in the partial derivatives of 
4>{%, y) will converge. Hence the conditions of Theorem V are satisfied. 

Conversely, if (35) converges, (36) will converge, and so z t will lie within 
or on the circle of convergence of f(z) with center a t . For the sum of m terms 
in (36) represents a root of the sum of corresponding terms in (35), and, if a 
matrix approaches a limit, we have shown that its roots must approach 
definite limits. Furthermore 

+ (<*<, *i) =/(*,) +/'(«<) (*<-«,) + • • • • +/"(«<) {Zi 7" i)m +■■■■ =/(*<)• 



m 



Also, if 






dx a dy e 
it is readily seen that 

^a.o(o«!*i) =/"(«<), ^«,/)(a*,«*)=0, /3=#=0. 
Hence (27) gives the same expression for <fy(A,Z) and for /(Z). Therefore 

f(Z) =f(A) +f(A) (Z-A) + .... +f m (A) ( Z ~ A ) m + (37) 

7/2' 

We may call this the expansion of f(Z) in the neighborhood of the matrix A. 
Theobbm VI. The Taylor's series (37) for f(Z) is valid for any matrix 
Z commutative with A if each root of Z lies within a circle with center at the 
corresponding root of A, in which f(z) is analytic. 

This theorem enables us to determine an analytic function of a matrix Z 
by means of a power series 

c 1 + c a (Z-A) + .... +c m (Z-A) m + ...., 

and its analytic continuations just as we determine an analytic function of the 
complex variable z by the series 

Ci + c 2 (z— a) + .... +c m (z— a) m + 
and its analytic continuations. The principal difference in the two cases is 
that Z and A must be commutative and that the region of convergence of the 
matrix series consists of m circles (one for each distinct root of A) instead of 
one. On a range of commutative matrices the function f(Z) is one valued if 
each root of Z is restricted to a region of the complex plane (not necessarily 
the same for different roots) in which the corresponding branch of f(z) is 
one valued. 



